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ABSTRACT 

We calculate light curves produced by r-modes with small azimuthal wavenumbers, 
m, propagating in the surface fluid ocean of rotating neutron stars. We include rela- 
tivistic effects due to rapid rotation, and propagate photons from the stellar surface 
to a distant observer using the Schwarzschild metric. The wave motions of the surface 
r-modes are confined to the equatorial region of the star, and the surface pattern of 
the temperature variation can be either symmetric (for even modes) or anti-symmetric 
(for odd modes) with respect to the equator. Since for the surface r-modes the oscil- 
lation frequency in the corotating frame of the star is much smaller than the rotation 
frequency, Q, we employ the approximation in which the oscillation frequency in the 
inertial frame, a, is given by a — —mQ. We find that the even, m — 1 r-mode pro- 
duces the largest light variations. The dominant Fourier component in the light curves 
of these modes is the fundamental having a = —Q, and the first harmonic component 
having er = — 2Q is always negligible in comparison. The dominant Fourier component 
of the even, m — 2 r-modes is the first harmonic. Although the odd r-modes produce 
smaller amplitude light variations compared with the even modes, the light curves of 
the former have a stronger first harmonic component. If both m = 1 and 2 r-modes 
are excited simultaneously, a rich variety of light curves is possible, including those 
having an appreciable first harmonic component. We show that the phase difference, 
5 — 5e, between the bolomctric light curve and that at a particular photon energy can 
possibly be used as a probe of the stellar compactness, R/M, where R and M are the 
radius and mass of the star. 
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1 INTRODUCTION 

It is now widely accepted that millisecond oscillations during thermonuclear X-ray bursts (hereafter, burst oscillations) on 
accreting, weakly magnetized neutron stars in low mass X-ray binary (LMXB) systems are produced by spin modulation of a 
slowly moving, non-uniform brightness pattern on the stellar surface (see Strohmayer et al. 1996; Strohmayer & Bildsten 2004, 
van der Klis 2004) . In addition to revealing the spin frequency these oscillations encode information about the global properties 
of the neutron star. A number of attempts have been made to model burst oscillations and thus to infer these properties (e.g., 
Nath, Strohmayer & Swank 2001; Muno, Ozel & Chakrabarty 2002; Poutanen & Gierliriski 2003, Bhattacharyya et al. 2005), 
however, uncertainties associated with the exact nature of the brightness pattern makes a unique interpretation of the data 
problematic. 

As of February, 2005, the number of LMXBs that exhibit burst oscillations amounts to 13, and two of these contain 
millisecond X-ray pulsars (see, e.g., Chakrabarty et al. 2003; Strohmayer et al. 2003). For most burst oscillations a single 
oscillation frequency is found that displays a positive frequency drift A^burst of a few Hz from the start to the tail of a burst. In 
the decaying tail the frequency often approaches an asymptotic value (e.g., Strohmayer & Markwardt 1999). The asymptotic 
frequencies are stable to a few parts in 1000 in bursts observed over several years from any given source (Strohmayer et al. 
1998b; Muno et al. 2002a). Background subtracted light curves for the burst oscillations are in most cases well fitted by a 
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single sinusoid, the frequency of which we will call the fundamental frequency. Harmonic content in the light curves is generally 
small (Strohmayer & Markwardt 1999; Muno et al. 2002b). Indeed, to date there exists only one strong detection of a first 
harmonic, in the millisecond X-ray pulsar XTE J1814-338 (Strohmayer et al. 2003, Bhattacharyya et al. 2005). On average 
the rms amplitudes of the burst oscillations are 2%-10%, and generally increase with photon energy (Muno et al. 2003). Muno 
et al.(2003) also explored the energy dependence of the oscillation phases of burst oscillations. They found either constant 
phase lags with energy or in some cases evidence for modest hard lags, behavior that is opposite to what one would expect 
from relativistic motion of a simple hot spot on the stellar surface (e.g., Ford 1999). Interestingly, where phase lags have been 
measured in the accreting millisecond pulsars they all show soft lags, consistent with a rotating hot spot model (Cui, Morgan 
& Titarchuk 1998; Poutanen & Gierlihski 2003). 

As noted above, a non-uniform brightness distribution on the surface of rotating neutron stars is almost certainly re- 
sponsible for the burst oscillations, and several models have so far been proposed to explain the inhomogeneous pattern. As 
a typical model, Strohmayer et al. (1997) suggested the oscillation originates from hot regions (hot spots) in a thermonuclear 
burning layer that expands and decouples from the neutron star at the start of the burst. Although a simple hot spot appears 
to be consistent with oscillations observed during the onset of some bursts (see Strohmayer et al. 1998a), further calculations 
suggested that a rigidly rotating, hydrostatically expanding burning layer produces too small a frequency drift (Cumming & 
Bildsten 2000, Cumming et al. 2002). 

There have been numerous efforts to calculate light curves produced by hot spots on the surface of rotating neutron stars 
(e.g., Pechenick et al. 1983, Chen & Shaham 1989, Strohmayer 1992, Page 1995, Braje et al. 2000, Weinberg et al. 2001, Muno 
et al. 2002, Nath et al. 2002, Poutanen & Gierliriski 2003, Bhattacharyya et al. 2005). In the modeling of the light curves many 
parameters are employed, for example, the number of hot spots, their location with respect to the rotation axis, their angular 
sizes, the inclination angle between the line of sight and the rotation axis, the angular dependence of the specific intensity, 
the temperature of the spots, the temperature contrast between the spots and their surroundings, the mass, radius, and 
spin frequency of the star. As noted above, burst oscillation light curves are remarkably sinusoidal. Although the light curve 
amplitudes are dependent on the relativistic parameter R/M, which is a parameter of particular importance, it is obvious 
that there remain parameters other than R/M that affect the amplitudes, such as the angular sizes and locations of the hot 
spots and the inclination angle (e.g., Pechenick et al. 1983, Weinberg et al. 2001, Muno et al. 2002b). In this sense it is quite 
important to find other observational quantities that can be used to obtain independent constraints on these parameters. 
The energy dependencies of phase lags are an example of potentially important observational quantities that carry useful 
information about the neutron star. Taking account of the special relativistic Doppler boosting of X-ray photons due to rapid 
rotation of the star (Ford 1999, Weinberg et al. 2001), the hot spot model can produce hard leads more or less consistent 
with those observed in the millisecond X-ray pulsar SAX J1808.4-3658 and ICR J00291+5934 (Poutanen & Gierlihski 2003; 
Galloway et al. 2005), but the model appears to be at odds with the suggestion of hard lags in some burst oscillations (Muno 
et al. 2003). Although it may be possible to produce hard lags via a hot Comptonizing corona, it is worth exploring in more 
detail whether more complex brightness patterns (such as oscillation mode patterns) on the stellar surface might be consistent 
with the production of hard lags. 

By inferring the hot spot size using the observationally obtained blackbody flux Ftt, Poutanen & Gierlihski (2003) have 
succeeded in estimating the radius of the millisecond X-ray pulsar SAX J1808. 4-3658 assuming an appropriate range of mass 
for the star. Recently, Strohmayer et al. (2003) have found the first harmonic component in the light curves of burst oscillations 
from the accreting millisecond X-ray pulsar XTE J1814-338, and Bhattacharyya et al. (2005) obtained a useful constraint on 
the ratio R/M assuming a single hot spot model. 

As alluded to above, the hot spot model is not always successful in producing the observed light curves of the burst 
oscillations. As discussed by Muno et al. (2002b), a small hot spot in a completely dark background is likely to produce light 
curves with too large an amplitude as well as having strong harmonic components not seen observationally. To suppress the 
large amplitudes and the harmonics one has to assume, for example, a hot spot of angular size a ~ 7r/2 or almost completely 
antipodal hot spots, the assumptions of which are not necessarily physically well motivated. It is therefore important to pursue 
alternatives to the hot spot model. 

In this paper, we calculate light curves of rotating neutron stars assuming the surface temperature perturbation pattern 
caused by surface r-modes propagating in the fluid ocean of a mass accreting neutron star (Lee 2004). It was Heyl (2004) who 
first argued that the r-modes with low m may play a possible role for burst oscillations. The surface r-modes are rotationally 
induced oscillations with very low frequencies in the corotating frame of the star. The r-modes have dominant vorticity 
(V x v') r and are basically just large scale vortexes (see Spitkovsky, Levin, & Ushomirsky 2002). If we assume axisymmetry 
of the equilibrium configuration of the rotating star, the time dependence of the temperature perturbation may be given by 
e lCT *, where a denotes the oscillation frequency in the inertial frame. If we let u> and fl denote the oscillation frequency in the 
corotating frame and the spin frequency of the star, respectively, we have a = to — mfl, which leads to a w — rriQ for the 
surface r-modes, since \ui\ « \mQ\, where m is an integer representing the azimuthal wave number around the rotation axis. 
The time dependence of the perturbation is thus to good approximation given by e _lmnt . Since only the oscillation modes 
having low m will produce light variations of appreciable amplitude, in this paper, we only consider the surface r-modes with 



Light curves of oscillating neutron stars 3 



m = 1 or m = 2. We briefly describe the method of calculation in §2. In §3 we discuss some of our numerical results with 
relevance to burst oscillations. §4 is for a discussion, and we give a summary of our principal conclusions in §5. 



2 METHOD OF SOLUTION 

The method of calculating light curves of rotating neutron stars is similar to those discussed by Pechenick et al. (1983), Page 
(1995), and Weinberg et al. (2001). Instead of presuming the presence of a hot spot on the surface of a rotating neutron star, 
we assume that the temperature varies accross the surface with the angular distribution appropriate for a surface r-mode 
propagating in the fluid ocean of the star. The surface r-modes, which have very low oscillation frequencies in the corotating 
frame of the star, can plausibly be excited during thermonuclear bursts (Lee 2004). 

To describe oscillations of rotating neutron stars, we employ an xyz coordinate system whose origin is at the center of 
the star with the rotation axis along the z direction. For convenience, we also assume that the observer is in the x-z plane. 
In this coordinate system, the oscillations of rotating neutron stars are described by employing spherical polar coordinates 
(r,9,<p) with the axis 9 = being the z-axis. Then, the displacement vector £(r,8,(f),t) is given by 

oo 

' ^SuYre^, (1) 



2=1 



*' = VA H h ia Y U + Tv --a^ ) e°\ (2) 



2 = 1 



d 



OC / 



and the Lagrangian perturbation of the temperature is given by 

OC 

8T = ^5T lj Y l ™e ia \ (4) 

2 = 1 

where a is the oscillation frequency in the inertial frame, and Y t m denotes the spherical harmonic function with degree I and 
azimuthal wave number m, and lj = \m\ + 2(j — 1) and lj = lj + 1 for even modes, and lj — m + 2(j — 1) + 1 and lj = lj — 1 for 
odd modes, and j = 1, 2, • • •. The functions £(r, 9, <f>, t) and ST are obtained as a solution to linear differential equations for 
oscillations of a rotating star, and we need to specify a normalization condition to determine the amplitudes. Since the surface 
r-modes have dominant toroidal components TJ/ , we will employ the amplitude normalization given by |Re[iT ; / (R)/R]\ = 1. 

3 1 

The surface temperature of a rotating star that pulsates in various oscillation modes may formally be given by 



T{R,0,<j>,t) =T + Rei ^C a 



I_2>1 

where To is the mean surface temperature, a indicates a combination of indices that distinguish various oscillation modes, 
and a a and C a denote, respectively, the oscillation frequency and a complex constant representing the amplitudes of the 
oscillation mode a. Note that by introducing C a in equation (5) we mean \Ke[iT" (R) / R]\ = \C a \. Using the oscillation 
frequency to = a + mQ. in the corotating frame, we have 

mcj) + a a t = m<t> + <*> a t, (6) 

where <j> = <j> — fit is the longitude in the corotating frame. 

To consider photon trajectories around a neutron star, it is convenient to introduce another coordinate system x'y'z' , 
for which the origin is at the center of the star, and the z' axis is pointing to the observer and y = y' . If we let i denote the 
inclination angle between the z axis and the z' axis, we have 

x = x cos i + z sin i, y = y' , z = — x sin i + z cos i, (7) 

and it is easy to derive the relation between the two spherical polar coordinate systems (r, 6, <fi) and (r', 8' , 4>'), where the axis 
of 9' = is the z axis. 

If we assume the Schwartzschild metric around the neutron star with mass M and radius R, the number flux, per photon 
energy Eoo, of thermal photons emitted from the stellar surface and received by a distant observer may be given by (e.g., 
Page 1995) 

1 d 2 iV(goo) _ 2tt Rl f2 [\ ^ 1 _ 2tt Rl El 

S dt^dE^ ~ (?h* °° J 9 27V efe-« E oo/k B T _ ! - C 2 h 3 D 2 ee -. EoD /k B T _ 1 <"S ee (t oo) , (tS) 
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and integrating equation (8) with respect to the variable Eoo leads to 

1 dN(E 
S dtoc 
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(9) 



where D is the distance to the observer, T is the surface temperature, k B is the Boltzmann constant, c is the velocity of light, 
and h is the Planck constant, and we have neglected, for simplicity, the interstellar absorption effect and the dependence of 
the effective area S of the detector on the photon energy. Here, 

-Roc = e~ n R (10) 
is the radius of the star seen by a distant observer, and dt x — e~ v dt denotes the increment of the coordinate time too, where 



e" = y/l ~ r g/R (11) 

and Rg = 2GM/c 2 is the Schwarzschild radius, and dt is the increment of the proper time at the surface of the star. Note 
that floe = e v Q since Qoodtoo = Sldt, where Qoo is the angular spin frequency of the star seen from infinity. If we let E e and 
Eoo denote, respectively, the photon energy in the corotating frame of the star at the surface and the photon energy recieved 
by a distant observer, we have 

E e = fe^Eoo, (12) 

where the factor /, representing the Doppler shift in the photon energy from e~ v Eoo in the non-rotating frame at the stellar 
surface to E e in the co-rotating frame, is given by 

/ = 7 (1 -vo/c), (13) 

where v = tt x R = RQsinOe,/, with e<f, being the unit vector in the azimuthal direction and o is the unit vector along the 
trajectory of a photon at the surface, and 7 = 1/ yl — v' 2 /c 2 . Note that the four vector of a photon in the non-rotating frame 
may be given by o M = e~ v E 00 (l, o). Using the angle C,{0') between the local normal and the photon trajectory, at the surface, 
reaching the observer at infinity, the variable q in equations (8) and (9) is defined by 

g^sinC, (14) 
and the relation between 0' and q is determined by the integration given by 

r R G /2R , 

m = / ==^= , —■ (15) 



A. - Ra/R) Rg/2R - (1 - 2u)itV 

The maximum value f^ax occurs when q = 1 and is shown as a function of the radius R in Figure 1 for M = IAMq. Since o 
at the stellar surface is given by 

o = sin(6*' — £) cos <f>'i' + sin(#' — £) sin <j> j' + cos(0' — Qk' , (16) 
where i' , j' , and k' are the orthonormal vectors parallel to the x' , y' , and z' axes, respectively, we have 

° - = sin(0' — () (cos <j) sin <f> — sin cos 4>' cos 1) — cos(0' — sin <fi sin i. (17) 

Note that without light bending so that 0' — we simply have v ■ o — —QRsin<j)smism0. The aberration effect may be 
given by 

COS ( e = — — COS C, (18) 

where ( e is the angle measured in the corotating frame of the star. 

For waves whose <j> and t dependence is given by e ™*'+ 1<Tt ) it is easy to incorporate the time-delay effect that arises because 
of the difference in the path lengths of photons emitted from different points on the surface of the star. Since the travel time 
of a photon from a point having the impact parameter b — qRoo at the surface of the star to an observer at infinity is given 
by 

= — r 1 (1 - (19) 

c Jr y/l-(l- R G /r){b/rY \ r J R G 
the time-lag may be defined by 

Af oo (&)EEf oo (6)-t oo (0). (20) 
Replacing too by toe — Atoo(fe), we rewrite the time and cj> dependence e lmi *+ lCT * Q f the perturbations as 

^im(f>-\-i<Tt ^ — ie 77 a At x (b) (21) 
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Table 1. Oscillation frequencies lo and a/m for the fundamental (n = 0) r-modes for the mass-accreting 
envelope model with M = O.lMEdd f° r M = 1.4Mq and R = 10km. Here the frequencies Co, a, and Q are 
normalized by 



k 


m 


n = 


0.2 

a/m 


Q = 

Co 


0.4 

a/m 





1 


5.31(-3) 


-0.1947 


5.31(-3) 


-0.3947 




2 


1.02(-2) 


-0.1949 


1.04 (-2) 


-0.3948 


1 


1 


3.24(-3) 


-0.1968 


3.22(-3) 


-0.3967 




2 


6.32(-3) 


-0.1968 


6.36(-3) 


-0.3968 



so that we can take account of the time-lag effects in light curve calculations (see e.g., Cadeau et al. 2005). 

For the light curves G(t) or Ge^ (t), we calculate the discrete Fourier transform cij (j = —N/2, ■ ■ ■ , N/2 — 1) defined by 

JV-l 

aj = G(t k )e 27,lf > tk (22) 

k=0 

where N is the total number of sampling points in the time span AT, and tk = kAT/N and fj — j/AT, and \aj\ = \ct-j\ for 
a real function G(t). To indicate the amplitude of the light variation with the frequency fj, we use the quantity defined by 

A(ojj) =2\a,j\/N, (23) 

where ujj = 2nfj . 



3 NUMERICAL RESULTS 

The surface r-modes propagating in the fluid ocean of a rotating neutron star may be classified in terms of three quantum 
numbers (m,k,n), where m is the azimuthal wave number around the rotation axis, k is the number of latitudinal nodes, 
and n is the number of radial nodes of the eigenfunctions. For the quantum number k, we count the latitudinal nodes of the 
eigenfunctions Re[ST(R, 8)], and for the quantum number n we count the radial nodes of the eigenfunction iTi^r). In this 
paper, we exclude the r-modes of I' — \m\ for light curve calculations, since their frequency changes during bursts are much 
larger than is suggested by observations (e.g., Cumming & Bildsten 2000, Lee 2004). Then, the r-modes with the quantum 
number k = 0, 1, 2, • • • correspond to the r-modes of I' = \m\ + 1, \m\ + 2, \m\ +3, • • •, respectively, and the corotating frame 
frequency ui decreases in this order for given m, n, and SI. The r-modes with even k are denoted "even" modes, and have 
eigenfunctions 5T(R, 6) that are symmetric about the equator of the star, while the r-modes with odd k are called "odd" 
modes and are antisymmetric about the equator. 

For light curve calculations in this paper, we consider only the fundamental r-modes with no radial nodes (n = 0) of 
the eigenfunction iTi j , because they are the most strongly excited in thermonuclear bursts (Lee 2004). In addition, since the 
r-modes associated with high m and k are unlikely to produce light variations with appreciable amplitudes, we restrict our 
numerical analysis to the fundamental modes with m = 1 or m = 2 and with k = or k = 1. Employing these simplifications, 
and assuming two r-modes with different m may be excited, we can express the surface temperature for light curve calculation 
as 

T(R, 6, t)=T + Re [dST ai (R, 6») e l W+ CT «i i ) + C 2 e x ST a2 (R, e)e l(2 ' l ' +a ^ t) ] , (24) 
where 

^(^^^(^(M), ST a *(R,e) = Y J 5T?;{R)Y?.{0,Q), (25) 
j>i j>i 

and Ci = C ai and Ci = C a2 are here real constants for the amplitudes of the modes, and \ is a real parameter, giving the 
phase difference between the two r-modes. If we define a = (m,k,n), we have ai = (1,0,0) and «2 = (2,0,0) for the even 
modes, and ai = (1, 1,0) and 02 = (2, 1,0) for the odd modes. For the oscillation frequencies, we adopt the approximation 
that a ai = — SI and a a2 = —2SI. The parameters we need to calculate the function G(t) defined by equation (9) are the 
mass M, the radius R, the angular rotation speed SI of the star, the oscillation amplitudes Ci and C2, the phase difference 
X between the two modes, and the inclination angle i. For the function Gs 00 (t) defined by equation (8), we need in addition 
the surface temperature To and the photon energy Eoo- 

For the r-mode calculation, we use a mass-accreting envelope model with the rate M = O.lMEdd for a neutron star having 
M = 1.4M© and R = 10 6 cm, where Msaa = 4arcR/K e with K e being the electron scattering opacity. This envelope model is 
a fully radiative model with no convective regions in it, and the detail of the envelope calculation is given in Strohmayer & 
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Lee (1996) and Lee (2004). To calculate the modes propagating in the mass-accreting envelope of a rotating neutron star, we 
follow the method described by Lee & Saio (1987, 1993). Note that the formulation employed for the envelope and oscillation 
mode calculations is Newtonian, and no general relativistic effects are included. In Table 1, the oscillation frequencies ui) and 
a/m arc tabulated for the fundamental r-modes with low m and k for f2 = 0.2 and Cl — 0.4, where the frequencies w, a, 
and Q. are normalized frequencies by \J GM/ R A . It is important to note that the corotating frame oscillation frequency uj is 
almost insensitive to Q for rapidly rotating stars (Lee 2004), and that for a given k the inertial frame pattern speed a/m 
for m — 2 is nearly equal to that for m = 1, which may justify the simplification of a = —mi~l. The temperature variations 
8T(R,9)/To caused by the surface r-modes for Q — 0.2 are shown in Figure 1, where panels (a) to (d) are for a = (1,0,0), 
(2, 0, 0), (1, 1, 0), and (2, 1, 0), respectively, and the solid and the dashed lines denote respectively the real and imaginary parts 
of the function 8T(R,8)/To. Here, the amplitude normalization is given by IRefiTJ^ (i?) / 7?] | = 1 at the surface of the star. We 
note that the functions 5T(R,8)/To and iTi 1 (R)/R have comparable amplitudes to each other at the stellar surface, and that 
the imaginary part of 8T(R,6) /To is much smaller than the real part. The 9 dependence of the function ST(R,9)/To for the 
even (odd) r-mode with m = 1 is quite similar to that for the even (odd) r mode with m — 2 except for the fact that the 
maxima of the functions for m = 2 are larger than those for m = 1. In this paper, the temperature perturbations ST(R,8) 
calculated for R — 10km will be used for different values of R for simplicity. Applying the approximation a — —mQ, we have 
m<f> + a a t — m<\> — mflt, and we call the Fourier amplitudes A(u>j) at ujj = f2 and uij = 2f2 the fundamental and the first 
harmonic components, respectively. 

Since we are applying a linear theory of stellar pulsations, we cannot self-consistently determine the amplitudes of the 
modes, which forces us to treat the oscillation amplitudes as parameters in the modeling. For simplicity, we use Ci = 0.2 for 
m = 1 and C 2 = 0.2 for m = 2 for the mode amplitudes throughout the following calculations. Since the toroidal component 
of the displacement vector is dominant for the r-modes, the velocity, v, of a surface fluid element may be approximately given 
by v ~ \ujiT h (R)\ = Ry/GM/R 3 \Cu\ x \iT h (R)/R\. Since |w| <D.01, we have the fluid velocity v <10~ 3 c for \iT h (R)/R\ ~ 0.2 

if we adopt the stellar parameters M = 1.4M Q and R = 10km, for which RW GMjR? ~ 0.5c. If we use the sound velocity 
v s = ksT/fiH at the surface, we have v s ~ 10~ 3 c for T = 10 7 K and fi = 0.5, where H and denote the hydrogen mass 
and the mean molecular weight, respectively. For the value of Ci = 0.2 (or C2 = 0.2) for the amplitude parameter, the fluid 
velocity at the surface is much smaller than the velocity of light and is comparable to or less than the local sound velocity at 
the surface. 



3.1 Light curves produced by a single r-mode with low m 

Since the imaginary part of ST(R,9)/Tq is much smaller than the real part, assuming that only a single r-mode of low m is 
excited, we may simply write the surface temperature as 



T = To [1 + C{0) cos(m0 - mfii)] , 



(26) 



where C(6) — Re[CjST(R,0)/To] and j = 1 or j = 2. Note that setting the origins of times t and £oo appropriately, we can 
assume Qt — Qactoo. The function G(t) is given by 



G(t) = G +Gi cos mm + G\ sin mQt + G c 2 cos 2mfit + G 3 2 sin 2mQ.t H , 

where 



Go 



2qdq 



64' 1 + 1.5C 2 ((9) 



2tt 



P 



2qdq 



1 + 



1 + 



C 2 {6) 



C 2 {6) 



cosm(4> + A(j>), 
sinm(0 + A0), 



J 2qdqJ ^-lf- 3 C 2 (e)cos2m^ + A4>), 



-f~ C {0) sin 2m(0 + A0), 



(27) 

(28) 
(29) 
(30) 

(31) 
(32) 



where A<p = f2oo At 00(b), and we have replaced too by too — Atoo(b) to calculate the light curves at too at the observer. It may 
be important to note that the coefficients Gf and G 2 vanish if we assume / = 1 and A<j> = 0. Equation (27) indicates that 
there appears a periodic component having the frequency 2vn£l unless G 2 and G 2 vanish simultaneously. The terms associated 
with the frequency mfl, for example, can be rewritten as 
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GI cos mQt + G{ sin mQt — G\ cos (mQt — 8) , (33) 

where Gi = ^(Gl) 2 + (GI) 2 , tan 5 = Gi/Gf, and the value of fit = (5 + 2ri7r)/m that maximizes the variation depends on 
the parameters R and i, where n is an integer. Since the imaginary part of ST(T, 6) is negligible compared with the real part, 
we have in a good approximation 

A(mQ) = Gi, A(2mQ) = G 2 , (34) 

where Gi = -J (Gj) 2 + (G|) 2 , and the Fourier amplitude -A(wj) is defined by equation (23). If we define the phase shift 8j, 
using the Fourier transform aj given in equation (22), as 

tanSj = Im(aj)/~Re(aj), (35) 

we have to good approximation 8 = 5j , selecting the integer j so that 2irfj = mQ. 

The quantities Gi and <5/(2rmr) defined by equation (33) are plotted as a function of the inclination angle i for the 
even (k = 0) r-mode with m = 1 in Figure 3, and for the even (k — 0) r-mode with m — 2 in Figure 4, where we assume 
Q = 0.2, and the dotted, short-dashed, solid, long-dashed, and dot-dashed lines indicate the radius R = 8, 9, 10, 15, and 
20km, respectively. The amplitude Gi vanishes at i = because the temperature variations are proportional to e 1 ™ 1 ^ where 
(j> = (j> for i = 0. For given i and R, the amplitude Gi for m = 1 is usually larger than that for m = 2. For a given inclination 
angle i, the amplitude A(mQ) = G\ decreases as the radius 7? decreases. This is because the maximum angle 8' msix increases 
with decreasing 7? so that almost the whole surface area of the star can be seen at any time by the distant observer. For the 
even modes, the amplitude Gi for a given R monotonically increases with increasing i. As shown by Figures 3 and 4 the phase 
shift 8/(2miv) for the even r-modes are negative. The absolute value \8/(2mn)\ for m = 1, which stays small, increases with 
decreasing R, and it is rather insensitive to the inclination angle i for a given 7?. The absolute value |<5/(2m7r)| for m = 2 
also increases as R decreases and becomes as large as |<5/(2m7r)| ~ 0.2 for the smallest R considered here. It shows a weak 
dependence on i, particularly for the radii R ~ 10km. 

The quantities Gi and 8/(2rmr) are given as a function of the inclination angle i for the odd (k — 1) r-mode with m = 1 
in Figure 5, and for the odd (k = 1) r-mode with m = 2 in Figure 6, where Q — 0.2, and the different curves have the same 
meaning as in Figure 3. For the odd modes, the amplitude Gi vanishes at i = 90° also, because the temperature variation 
pattern is antisymmetric about the equator of the star. As in the case of the even r-modes, the amplitude Gi for the odd 
r-modes decreases with decreasing R for a given i. Although 8/(2rmr) for larger radii i?^>15km stays constant with varing i, 
it shows a strong dependence on i for smaller radii R <10km. 

Figure 7 shows the amplitude ratio G2/G1 for the even r-mode with m = 1 (upper panel) and for the odd r-mode with 
m = 1 (lower panel) as a function of i. We find the ratio is at most of order ~0.01 for the even m = 1 r-modes, which 
indicates that the first harmonic component cannot be significant in the light curves produced by a single even r-mode of low 
m. For the odd m = 1 r mode, on the other hand, the first harmonic component has appreciable amplitude compared to the 
fundamental, particularly when i ~ 90°. Actually, the ratio G2/G1 diverges as i increases to 90°. 

As examples, we show as a function of Q.t/2-K the light curves produced by the even, m = 1 r-mode in Figure 8, and by 
the odd, m = 1 r-mode in Figure 9, where we assume f2 = 0.2 and i = 60°. Figure 10 shows the light curves generated by the 
odd, m = l r-mode as a function of Qt/2n for Cl — 0.2 and R = 10km, where the dotted, short-dashed, solid, long-dashed, 
and dot-dashed lines indicate the inclination angle i — 10°, 30°, 50°, 70°, and 90°, respectively. At i = 90°, only the first 
harmonic component associated with 2mf2 appears. 

3.2 Dependence on f2 

As Q, increases, the amplitudes of the r-modes are more strongly confined to the equatorial region of the star, which makes the 
amplitude Gi larger for the even r-modes but smaller for the odd modes. Note that the frequency uj is only weakly dependent 
on Q, for the modes. To see the effects of rapid rotation, we calculate Gi and 8 for the even, m = 1 r-mode for Cl — 0.4, and we 
plot the results in Figure 11, where the different curves represent different radii as in Figure 3. We note that the magnitude 
of the phase shift 8 is larger for Cl — 0.4 than for Q, — 0.2 although 8 is still only weakly dependent on the inclination angle i 
for the mode. 

3.3 Light curves produced by two r-modes 

Although the even, m = 1 r-mode will produce the largest light variations of the modes considered here, the mode is not 
likely to produce a first harmonic component A(2Q) of appreciable strength. In order to account for observed light curves that 
contain a substantial first harmonic component, we may need to assume that the r-modes with m = 1 and m = 2 are excited 
simultaneously. As an example, we calculate the light curves G(t) produced by the two simultaneously excited even r-modes 
with m = 1 and m = 2, assuming Gi = C2 = 0.2, x — 0, i = 90°, and Cl — 0.2. We plot G(t) as a function of Qt/2n in 
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Figure 12 and the Fourier amplitudes A(Q) and A(2£l) as a function of the radius R in Figure 13, where the different curves 
in Figure 12 have the same meaning as in Figure 3. Obviously, using two r-modes with different m's, we can produce a variety 
of light curves, including those with a harmonic content similar to that seen from XTE J1814-338 (Strohmayer et al. 2003). 
It is important to note that, even if we abandon the simplification given by a ai = — and a a2 = — 20,, we have in a good 
approximation 2a ai = a a2 as shown by Table 1. 



3.4 Phase Lags 

The amplitudes of the light curves produced by the surface r-modes are inevitably related to the parameters d's as well as 
to the ratio R/M. Since the oscillation amplitudes are difficult to determine within the framework of a linear theory of stellar 
pulsations and the ratio R/M itself is one of our main parameters to determine observationally, it would be useful to find an 
additional observable that can be used as an indicator for the ratio R/M. Considering the even r-mode with m = 1, which 
is most likely to be responsible for the light variations, we note that the phase shift 8, which is almost independent of the 
oscillation amplitude parameter Ci , is rather insensitive to the inclination angle i but is almost a monotonic function of the 
radius R (i.e., the ratio R/M) for a given i. If it is true that the observed light variations during burst oscillations are caused 
by an even, m = 1 r-mode, it would be possible to use the phase shift 8 as an indicator for the ratio R/M. To derive an 
observable using the phase shift 8, we calculate the function Gs^t) defined in equation (8) for appropriate values of To and 
E x , and Fourier analyze the function G_B 00 (t) to obtain the phase shift 5j at coj = O using equation (35), which we denote 
5e- In Figure 14, we plot as a function of the inclination angle i the difference 8 — 8e for the even, m = 1 r-mode (upper 
panel) and for the odd, m = 1 r-mode (lower panel) where we have assumed ft = 0.2, E x — E = lkev, fcsTo = 2.3kev, and 
the various curves represent different radii as in Figure 3. Since the difference 8 — 8e for the even r-mode with m = 1 is 
almost insensitive to i and its magnitude increases almost monotonically with the radius R for a given i, the difference may 
be useful as an observational indicator for the ratio R/M. On the other hand, the difference for the odd mode shows a strong 
dependence on i for the radius R <10km. Since the phase shifts 8 found in this paper are due mainly to the Doppler effects 
associated with rapid rotation of the star, the magnitude of the difference will be large for higher O for given M and R. 

As shown in Figure 14, the phase difference, 8 — 8e , for the even mode is always negative, but that for the odd mode 
can be positive for radii R <10km. Here, negative and positive values of 8 — 8e indicate hard leads and lags, respectively, 
for Eoc — Eo = lkev and UbTq = 2.3kev. To compare with observations more directly, it is convenient to give the phase shift 
difference, 8e — 8e , between the phase for Eoo — E and that for a particular E x = Eo. Examples are given as a function 
of the photon energy E in Figure 15 for the even r-mode with m = 1 for i = 60° and in Figure 16 for the odd r-mode with 
m = 1 for i — 30° (upper panel) and for i = 60° (lower panel), where we have used Eo = lkev and Qoo/2n — 400Hz. For 
the even mode we always have negative Se — 8e , indicating hard leads (soft lags), as suggested by Figure 14. On the other 
hand, for the odd m = 1 r-mode, we have both hard leads (soft lags) and hard lags, the latter of which occur at smaller radii 
(i.e., smaller R/M). Observationally, there is an indication for hard lags during some burst oscillations, however, others are 
consistent with no phase variations with energy (Muno et al. 2003). Since the even, m = 1 r-mode generates hard leads (as 
the hot spot model does), our results suggest that a single even r-mode is probably inconsistent with burst oscillations which 
show hard lags. However, it might be possible for an oscillation produced by contributions from both even and odd modes 
to match the observed phase lag properties of burst oscillations. We will explore detailed modeling and fitting, including the 
effects of the detection process, in a subsequent paper. 



4 DISCUSSION 

The surface r-modes considered in this paper belong to a subclass of the equatorial waves (Pedlosky 1987). In the equatorial 
f3 plane approximation, we employ a Cartesian coordinate system in which the coordinate origin is at the equator and the 
x, y, and z axes are in the eastward, northward, and upward directions, respectively. The wave equations in the /3 plane 
approximation are solved by applying the method of separation of variables, in which a separation constant A is introduced 
between (x, y, t) and z. The separation constant A is an eigenvalue that is determined by solving the wave equation in 
the vertical direction with appropriate boundary conditions, and reflects the vertical structure of the thin fluid envelope. 
Assuming the perturbed velocity v' y (x,y, z,t) is given by v' y oc e t ^ kx + ut ')^ (y)V (z) , we obtain for the equatorial waves the 
dispersion relation given by (Pedlosky 1987) 

\Cb 2 + k/Cj -k" 2 = (2j + l)V\, (36) 
where 

Cj = uj/((3 L e ), k = L e k, L e = y/NoD/fh, /3 = 2Q,/R, (37) 

and iVo is a characteristic value of the Brunt- Vaisala frequency, D is the depth of the fluid ocean, and the function ^(y) is 
given by 



Light curves of oscillating neutron stars 9 



tt(y) = Mr) = e-^/'fTjW/^i^, (38) 

where Hj(r) is the Hermite polynomial, r = A 1//4 t//L e , and j is an integer. Here, the function V(z) is determined by the wave 
equation in the vertical direction. The pressure perturbation for the equatorial waves for j > 1 is then given by (Pedlosky 
1987) 



V2^ d-h/y/X ^ 2 J 



(40) 



u-k/\f\ V 2 / [i + fc/v^. 
where is an arbitrary constant. Since the low frequency equatorial waves at rapid rotation have (Lee 2004) 
mNoD/R 
(2j + 1)VX 

where we have assumed k ~ m/i?, p'(a;, y, z, t) reduces to for j — 1 

„ _ Re {j^^ (r> + i) »V(„} , (41, 

and for j = 2 

?'(*,„,*,«) ~ -Re I ^T77TT7^ , V e^V^'V^l ■ (42) 



Knowing that the perturbed potential temperature ©' is approximately given by ©' = dp' /dz (Pedlosky 1987), we find that 
the latitudinal (y) dependence of <d'(x,y,z,t) well reproduces that of ST(R,6)/To in Figure 2 both for the even r-modes 
(j — 1) and for the odd r-modes (j = 2). Note that the solutions having the frequency (40) correspond to the type 2 solutions 
discussed by Longuet-Higgins (1968). 

The oscillation energy Ew of a mode observed in the corotating frame may be given by (e.g., Unno et al. 1989) 

f f ^*j»,f _iGMAM a 

E w = — £-£dM r =o/ , (43) 

2 Jam. k 

where 

AM. = \ I ^fdM r , (44) 
z Jam. rt 

and is the complex conjugate of £, and AM e is the envelope mass. If we let AEt and AE P denote the energies released during 
a burst and a persistent phase, respectively, we may assume a = AE p /AEb ~ 10 2 and AE P ~ eGMAM a /R, where e ~ 0.1 
is a factor representing the efficiency of energy transformation from gravitational to radiation energies and AM a = MAt p 
with M ~ 10 16 gs _1 being a typical mass accretion rate in close binary systems is the amount of mass accreting during the 
persistent phase At p ~ 10 4 s (see, e.g., Frank, King, & Raine 2002). We then have 

E w = ae~ 1 Q 2 {AM / AM a )AEb, (45) 
which may be rewritten with appropriate normalization as 

£w ~ 2 X 100 loij ^5xlO-«M s-iJ vlC^sJ WJ vOlJ 10-i0M o Ai5i " (46) 
where we have defined 

Li 

i? 2 



{AM ) = \f y-dM r , (17) 

Z JAAf c 

assuming £ = C| with £ being the displacement vector normalized as IRe^T;/ (R)/R]\ = 1. For the r-modes discussed in this 
paper, we have the ratio (AM ) /1O~ 1O M ~ 0.1 for the even r-modes, and (AM a ) /10~ 10 M Q ~ 0.01 for the odd r-modes, 
for f2 = 0.2 and AM e — 10 _1 Mq. Although the estimation given above contains many parameters, all of which are not 
necessarily well determined, we may expect that the amount of energy required to excite the modes to appreciable amplitudes 
is a modest fraction of that released in a burst. It is interesting to note that, although the odd r-modes need to have larger 
amplitudes to produce appreciable light variations than the even r-modes, the former are easier to excite than the latter for 
a given amount of energy released in a burst. 



5 CONCLUSIONS 

We have calculated light curves produced by the surface r-modes of a rotating neutron star, taking account of the effects 
of gravitational light bending, gravitational red shift, and the difference in the arrival times of photons traveling in the 
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static Schwarzschild spacetime, as well as the Doppler shift of photon energy due to rapid spin of the star. We find that the 
fundamental even r-mode with m = 1 produces the largest light variations. The light curves produced by a single even r-mode 
of a given m are dominated by the fundamental component with frequency mfi, although those produced by a single odd 
r-mode also contain a first harmonic component of appreciable amplitude as well as the fundamental component, the relative 
strengths depending on the inclination angle i. The phase shift 5 — 8e produced by the even m = 1 r-mode is only weakly 
dependent on the inclination angle i and is almost a monotonic function of R for a given i. The phase shift produced by the 
even r-mode with m = 1 is a hard lead. The phase shift produced by the odd, m = 1 r-mode, on the other hand, depends 
both on i and on R, and there exists a parameter space of R and i which produces hard lags. 

It is useful to translate the phase shift difference into a hard-lag measure, AiLag, which is given for 0.^/2% = 400Hz by 

2-7T 5e — 5e _ _ Se — Se ,, o1 

At ^ = ^^r = 2 - 5x ^r ms ' (48) 

where negative Afcag means hard leads. As indicated by Figure 15, to fit the observed hard leads ~ — 200/xs at Eoo ~ lOkev 
for SAX J1808. 4-3658 (Cui et al. 1998) in terms of the even m = 1 r-mode we need a little bit smaller values of the ratio R/M. 
In fact, Weinberg et al. (2001) have assumed M = 2.2M© and R = 10km for their fit, which is equivalent to R/Rg = 1.54, 
the value of which is substantially smaller than the smallest value R/Ra = 1-94 for M = 1.4M© and R = 8km we use in 
this paper. As suggested by Figure 16, however, the phase shift difference produced by the odd m = 1 r-mode may give a 
consistent value ~ — 200ps at E x ~ lOkev for M — 1.4M©, R — 10km, and i — 60°, although the strong saturation of AtLag 
above E >10kev is not necessarily well reproduced. 

Compared with the hot spot model, the wave model generally produces smaller amplitude light variations. This is mainly 
because the entire surface of the neutron star produces X-ray emission. On the other hand, purely sinusoidal light curves 
observed for most of the burst oscillators (Muno et al. 2002b) can be explained by assuming the dominance of a single r-mode 
of a given m, and the existence of an appreciable first harmonic component may be attributable to the coexistence of a mode 
with 2m, or perhaps a superposition of even and odd modes of different amplitude. As in the case of the hot spot model, the 
wave model also gives hard leads in most of the parameter space, which may be in conflict with the suggested hard lags seen 
in some of the burst oscillators, but might be applicable to the hard leads observed from the millisecond X-ray pulsars SAX 
J1808.4-3658, and ICR J00921+5934 (e.g., Cui et al. 1998; Galloway et al. 2005). The parameters employed in the light curve 
calculation include the mass M, the radius R, and the spin frequency Q of the star, the inclination angle i, the mode indices 
(m, k, n), and the oscillation amplitudes Cj. Even if we assume that only the fundamental n — r-modes with k — or k — 1 
play a role in producing the observed brightness variations and that the frequency that appears in the light curves is the 
spin frequency of the star, we still have several free parameters that affect the amplitude determination of the light curves. 
Perhaps the most crucial parameters are the oscillation amplitudes, Cj, themselves, which are difficult to determine within 
the framework of linear pulsation theory. Further theoretical studies combined with model fitting to the observed properties 
of burst oscillations will be definitely desired for a better understanding of the burst oscillations and the underlying neutron 
stars. 
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Figure 1. 9' mwc as a function of the radius R for the star with the mass M = IAMq 

REFERENCES 

Bhattacharyya, S., Strohmayer, T.E., Miller, M.C., & Markwardt, C.B. 2005, ApJ, 619, 483 
Braje, T.M., Romani, R.W., & Rauch, K.P. 2000, ApJ, 531, 447 
Cadeau, C, Leahy, D.A., k, Morsink, S.M. 2005, ApJ, 618, 451 

Chakrabarty, D., Morgan E.H., Muno, M.P., Galloway, D.K., Wijnands, R., van der Klis, M., & Markwardt, C.B. 2003, Nature, 424, 42 

Chen, K., & Shaham, J. 1989, ApJ, 339, 279 

Cui, W., Morgan, E.H., & Titarchuk, L.G. 1998, ApJ, 504, L27 

Cumming, A., & Bildsten, L. 2000, ApJ, 544, 453 

Cumming, A., Morsink, S.M., Bildsten, L., Friedman, J.L., & Holz, D.E. 2002, ApJ, 564, 343 
Ford, E.C. 1999, ApJ, 519, L73 

Frank, J., King, A., & Raine, D. 2002, Accretion Power in Astrophysics (3rd ed.; Cambridge: University of Cambridge Press) 

Galloway, D.K., Markwardt, C.B., Morgan, E.H., Chakrabarty, D., & Strohmayer, T.E. 2005, astro-ph/0501064| 

Gierliriski, M., Done, C, & Barret, D. 2002, MNRAS, 331, 141 

Heyl, J.S. 2004, ApJ, 600, 939 

Lee, U. 2004, ApJ, 600, 914 

Lee, U., & Saio, H. 1987, MNRAS, 225, 643 

Lee, U., & Saio, H. 1993, MNRAS, 261, 415 

Longuet-Higgins, M.S. 1968, Phil. Trans. R. Soc, 262, 511 

Miller, M.C., & Lamb, F.K. 1998, ApJ, 499, L37 

Muno, M.P., Chakrabarty, D., & Galloway, D.K. 2002a, ApJ, 580, 1048 
Muno, M.P., Ozel, F., k, Chakrabarty, D. 2002b, ApJ, 581, 550 
Muno, M.P., Ozel, F., & Chakrabarty, D. 2003, ApJ, 595, 1066 
Nath, T.R., Strohmayer, T.E., & Swank, J.H. 2002, ApJ, 564, 353 
Page, D. 1995, ApJ, 442, 273 

Pechenick, K.R., Ftaclas, C, & Cohen, J.M. 1983, ApJ, 274, 846 
Pedlosky, J. 1987, Geophysical Fluid Dynamics (2d ed.; New York: Springer) 
Poutanen, J., & Gierliriski, M. 2003, MNRAS, 343, 1301 
Spitkovsky, A., Levin, Y., & Ushomirsky, G. 2002, ApJ, 566, 1018 
Strohmayer, T.E. 1992, ApJ, 388, 138 

Strohmayer, T.E., & Bildsten, L. 2004, in Compact Stellar X-Ray Sources, ed. W.H.G. Lewin &: M. van der Klis (Cambridge: Cambridge 
Univ. Press) 

Strohmayer, T.E., Jahoda, K., Giles, B.A., & Lee, U. 1997, ApJ, 486, 355 
Strohmayer, T.E., & Lee, U. 1996, ApJ, 467, 773 
Strohmayer, T.E., & Markwardt, C.B. 1999, ApJ, 516, L81 

Strohmayer, T.E., Markwardt, C.B., Swank, J.H., k, in'T Zand, J. 2003, ApJ, 596, L67 
Strohmayer, T.E., Zhang, W., Swank, J.H., White, N.E., & Lapidus, I. 1998a, ApJ, 498, L135 
Strohmayer, T.E., Zhang, W., Swank, J.H., & Lapidus, I. 1998b, ApJ, 503, L147 

Strohmayer, T.E., Zhang, W., Swank, J.H., Smale, A., Titarchuk, L., Day, C, & Lee, U. 1996, ApJ, 469, L9 



12 U. Lee & T.E. Strohmayer 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

COS0 COS0 

Figure 2. ST(R, 0)/To as a function of cos 8 for Q = 0.2. Panel (a) for the fundamental even r-mode with m = 1, (b) for the fundamental 
even r-mode with m = 2, (c) for the fundamental odd r-mode with m = 1, (d) for the fundamental odd r mode with m = 2, where 
the solid and the dotted lines indicate the real and the imaginary parts, respectively. Here, the amplitude normalization for the mode is 
given by \Ke[iTi 1 (R) / R] \ = 1 at the surface. 
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Figure 3. Amplitude G\ and phase shift <5 as a function of the inclination angle i for the fundamental even r-mode with m = 1 (Ci 
C2 = 0, and x = 0), where Q = 0.2, and the dotted, short-dashed, solid, long-dashed, and dash-dotted lines are for the radii _R 
10, 15, and 20km, respectively. Here, the amplitude normalization for the mode is given by |Re[iTj (R)/R] \ = 1 at the surface. 




Figure 4. Same as Figure 3 but for the fundamental even r-mode with m = 2 (Ci = 0, C'2 = 



.2, and \ = 0). 




Figure 5. Same as Figure 3 but for the fundamental odd r-mode with m = 1 (Ci = 0.2, C2 = 



0, and X = 0). 




Figure 6. Same as Figure 3 but for the fundamental odd r-mode with m = 2 (Ci = 0, Ci = 0. 



2, and \ = 0). 




Figure 7. Ratio G2/G1 as a function of the inclination angle i for the fundamental even r-mode with m = 1 (upper panel) and for the 
fundamental odd r-mode with m = 1 (lower panel), where Q = 0.2, and the various curves represent different radii as in Figure 3. Here, 
the amplitude normalization for the mode is given by \Ke[iTi 1 (R) / R] \ = 1 at the surface. 
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Figure 8. G(i) as a function of Qt/2n for the fundamental even r-mode with m = 1 for i = 60°, where C\ = 0.2, C2 = 0, \ = 0, 
and Q = 0.2, and the various curves represent different radii as in Figure 3. Here, the amplitude normalization for the mode is given by 
\Re[iTi 1 (R)/R]\ = 1 at the surface. 
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Figure 9. Same as Figure 8 but for the fundamental odd r-mode with m = 1. 
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Figure 10. G(t) as a function of Qt/2n for the fundamental odd r-mode with in = 1 for R = 10km, where Ci = 0.2, C2 = 0, \ = 0, and 
Q = 0.2, and the dotted, short-dashed, solid, long-dashed, and dash-dotted lines indicate the inclination angle i = 10°, 30°, 50°, 70°, 
and 90°, respectively. Here, the amplitude normalization for the mode is given by \Ke[iTi 1 (R) / R}\ = 1 at the surface. 
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Figure 11. Same as Figure 3 but for Q = 0.4. 
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Figure 12. G(t) as a function of Qt/2ir for the fundamental even r-modes with m = 1 and m = 2 for i = 90°, where C\ = C'2 = 0.2, 
X = 0, and Q = 0.2, and the various curves represent different radii as in Figure 3. Here, the amplitude normalization for the mode is 
given by IRe^Tjj (R)/R]\ = 1 at the surface. 
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Figure 13. Fourier amplitudes A(Q) and A(2Q) as a function of the radius R for the fundamental even r-modes with m = 1 and m = 2 
for i = 90°, where C\ = C'2 = 0.2, \ = 0, and Q = 0.2, and the solid and dotted lines indicate A(Q) and A(2Q), respectively. Here, the 
amplitude normalization for the modes is given by |Re[iT; 1 (.R)/.R]| = 1 at the surface. 
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Figure 14. Phase shift difference (S — 5g ) / (2mir) as a function of the inclination angle i for the fundamental even m = 1 r-mode (upper 
panel) and for the fundamental odd m = 1 r-mode (lower panel), where Q = 0.2, Ci = 0.2, C2 = 0, and x = 0, and the various curves 
represent different radii as in Figure 3. Here, we have assumed kgTo = 2.3kev and Eoo = Eq = lkev to calculate the function Gj^ (t) 
and the phase shift 8e - Negative values of (<5 — 8 e q ) / (2mir) indicate hard leads. 
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Figure 15. Phase shift difference (8e — $E )/(2m7r) as a function of the photon energy E^o = E for the fundamental even m = 1 r-mode, 
where Q.oo/2-n = 400Hz, Ci = 0.2, C2 = 0, \ = 0, and i = 60°, and the various curves represent different radii as in Figure 3. Here, we 
have assumed kgTo = 2.3kev and E^ = Eo = lkev to calculate the function Gg^ft) and the phase shift difference (<5g — 8 e ) / (2rmr) , 
the negative values of which indicate hard leads. 
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Figure 16. Phase shift difference (Se — &E )/ (2rrt7r) as a function of the photon energy Eoo = E for the fundamental odd m = 1 r-modc 
for i = 30° (upper panel) and for i = 60° (lower panel), where f2oo/27r = 400Hz, Ci = 0.2, C*2 = 0, and \ = 0> anc i * ne various curves 
represent different radii as in Figure 3. Here, we have assumed kgTo = 2.3kev and E^, = Eq = lkev to calculate the function Gb^ (t) 
and the phase shift difference (5g — 5^ )/(2m7r), the negative values of which indicate hard leads. 



